Magnetoresistance and Magnetic Ordering Fingerprints in Hydrogenated Graphene 
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Spin-dependent features in tlie conductivity of grapliene, chemically modified by a random dis- 
tribution of hydrogen adatoms, are explored theoretically. The spin effects are taken into account 
using a mean-field self-consistent Hubbard model derived from first-principles calculations. A Kubo- 
Greenwood transport methodology is used to compute the spin-dependent transport fingerprints of 
weakly hydrogenated graphene-based systems with realistic sizes. Conductivity responses are ob- 
tained for paramagnetic, antiferromagnetic, or ferromagnetic macroscopic states, constructed from 
the mean-field solutions obtained for small graphene supercells. Magnetoresistance signals up to 
~ 7% are calculated for hydrogen densities around 0.25%. These theoretical results could serve as 
guidance for experimental observation of induced magnetism in graphene. 
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Introduction- Close to the charge neutrality point, 
two-dimensional graphene exhibits fascinating transport 
properties deriving from the massless Dirac fermions 
physics [J, These electronic properties can be fur- 
ther tuned and diversified through chemical functional- 
ization 0] or irradiation (defects formation) Indeed, 
strong modifications of sp^-bonded carbon network by 
grafting sp^-type defects is, for instance, an interesting 
strategy to induce an insulating behavior in graphene. 
Following this idea, recent attempts based on hydrogena- 
tion and fluorination were able to turn graphene into a 
wide band-gap insulator 0, Q . 

While ideal bidimensional graphene is nonmagnetic, its 
unidimensional derivatives (graphene nanoribbons) could 
exhibit magnetic ordered states due to the presence of 
edges 01 ■ Additionally, monovacancies or sp'^-type de- 
fects (such as those created by adsorbed atomic H) can 
introduce local sublattice imbalances in graphene, thus 
inducing magnetic moments 0, Q according to Lieb's 
theorem [lo| . Experimentally, the existence of intrinsic 
magnetism in graphene and graphite-based materials is 
strongly debated [III , [l^ . Some experimental works as- 
sume that the observed magnetic ordering is coming from 
graphene edges or grain boundaries 12|, a scenario that 
has been recently questioned by Sepioni and coworkers 
who found no magnetism even for high concentration of 
edge defects [l^l ■ Differently, in a recent experiment 14 1 , 
signatures of room temperature ferromagnetism in metal 
free untreated graphite were tentatively related to the 
presence of hydrogen impurities. We note that in all 
those experiments, the nature and quantity of defects re- 



main poorly characterized and mostly uncontrolled. It 
would be consequently highly desirable to clarify the 
conditions for observing spin-dependent transport finger- 
prints in hydrogenated graphene. 

In this Letter, spin-dependent transport properties of 
hydrogenated graphene-based systems with realistic sizes 
are explored theoretically. Firstly, a spin-dependent one- 
orbital Hubbard Hamiltonian, describing adsorbed hy- 
drogen impurities on a graphene sheet, is derived using 
first-principles calculations and solved at a mean-field 
level. The local spin texture around hydrogen adatoms 
is found to be strongly dependent on both the concen- 
tration and the specific adsorption sublattice of neigh- 
boring adatoms as previously reported [3l ■ The self- 
consistent Hubbard Hamiltonian is then implemented 
into a real space (Kubo)-transport method. Our theo- 
retical results reveal the suitable hydrogenation coverage 
of graphene to maximize spin-dependent conductivities, 
which translate into measurable magnetoresistance sig- 
nals in the diffusive regime. 

One-orbital mean-field Hubbard approximation- The 
one-orbital tight-binding description of graphene restricts 
to Pz orbitals centered on each carbon atom. The spin- 
related physics induced by Coulomb interaction is intro- 
duced in the calculation by means of the Hubbard model 
in its mean-field approximation 

H = i ^ 4<TCj,<T + f/5Z"'.t('^u) +"i,i("i,t): (1) 
<i,j>,o- i 

where t is the first-neighbors hopping term, (cj,cr) is 
the creation (annihilation) operator in the lattice site i 
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(j) with spin cr, U is the on-site Coulomb repulsion, and 
'T-i.ij 'T-i.t ^'^s the self-consistent occupation numbers for 
spin-down and spin-up electrons, respectively. The ratio 
U /t has been chosen to accurately reproduce various 
physical magnitudes extracted from first-principles 
calculations (as explained below). When a hydrogen 




FIG. 1. (color online) (a) Schematic representation of the sp^ 
hybridization breakdown upon hydrogen adsorption. In the 
present one-orbital tight-binding model, such a defect is equiv- 
alent to the creation of a monovacancy. (b) Single-vacancy 
supercell used in the mean-field Hubbard model calculations 
with periodic boundary conditions, (c) Supercell with two 
vacancies. The distance between two hydrogen defects is la- 
beled by d. The green colored region give the space location 
for the defects, to minimize the correlations between vacancies 
in neighboring cells. 



scopic ground state is given by the excess of magnetic 
moments on one specific sublattice, although S" = is 
the most likely value on simple statistical grounds (equal 
H occupation of both sublattices). In the following, two 
assumptions related to the macroscopic magnetic states 
are considered. The first one is the simplest: a param- 
agnetic (PM) state, which is a good approximation for 
sufficiently low H concentrations, or at sufficiently high 
temperatures, in which there is no correlation between 
the semi-local moments. This state is thus constructed 
from independent large supcrcells containing a single 
vacancy as shown in Fig IT] (b). The corresponding spin 
densities on each cell site is represented in FigI21[a,d) as 
a function of the distance to the center of the supercell 
for two different concentrations of defects. The result 
favorably compare with first principles calculations when 
U /t is chosen appropriately [l8| . 
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(H) atom is adsorbed on top of a carbon (C) atom, 
the sp^-symmetry is locally broken, and the electron 
from the C Pz orbital is removed from the tt bands to 
form a a bond with the H atom. Consequently, within 
the tight-binding approximation, such an adsorption is 
modeled by removing the corresponding electron and 
lattice site [as illustrated in FigH] (a)]. Figure [TJb) 
presents one of the supercells {N = 791 atomic sites) 
used in periodic boundary conditions calculations. This 
hydrogen defect induces zero-energy electronic states 
which are mainly localized around the impurity 
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When the Coulomb repulsion is switched on in the calcu- 
lation, these zero-energy states spin-polarize, leading to 
semi-local 5 = 1/2 magnetic moments with a staggered 
spin density mostly located on one sublattice d, 15 1 
[see Figl5] (a,d)]. When a finite concentration of de- 
fects is present, the magnetic ordering between these 
magnetic moments is dictated by the sublattices where 
they are located, being co-polarized or ferromagnetic 
(FM) for the same sublattice and counter-polarized or 
antiferromagnetic (AF) otherwise [Figl31^b,e)]. 

Magnetic states.- The total spin S of the macro- 



FIG. 2. (color online) Local (site) spin (s) versus r (the dis- 
tance to the center of the supercell) for various vacancy con- 
centrations: (a,d) rix ~ 0.1% and 0.3% in the paramagnetic 
case, (b,e) = 0.2% {d = 27.7A) and 0.4% (d = 18. lA) in 
an antiferromagnetic state, and (c,f) in a ferromagnetic ex- 
cited state within the mean-field Hubbard approximation [t 
(red) and J, (blue) electrons]. In all cases, the center of mass 
of the vacancies is at r = 0. 



The second assumption is more complex since pair- 
wise magnetic ordered states are investigated. The latter 
are constructed using supercells containing two vacan- 
cies on different sublattices [see FigH] (c)] far from the 
edges. When selecting one site on each sub-lattice [see 
Figlljb,e)], the magnetic moments couple antiferromag- 
netically and the macroscopic state generated by merging 
several supercells is an antiferromagnet (AF) with zero 
total spin (5* = 0). We explore the effect of the distance 
(d) between vacancies while varying accordingly the su- 
percell size. By changing simultaneously the distance and 
size, we always allow for a meaningful definition for the 
concentration based on only two vacancies per cell. 
One notes that for rix = 0.6% {d = 8.5A) the vacancies 
become very close to each other, thus allowing overlap- 



ping between magnetic moments and suppression of the 
local antiferromagnetic state We note that, when the 
vacancies are in the same sublattice, the ground state is 
ferromagnetic (FM), but we will not consider this case in 
this work. 

The AF ordering can be further tuned by applying a 
sufficiently large external magnetic field. This has been 
estimated to lie in the range B w 1 — 10 T for hydro- 



genated graphene nanoribbons 15j and similar distances 
between H adatoms. Since the magnitude of the local 
magnetic moments is lower in the case of graphene due 
to the absence of confinement, we expect the number to 
be smaller in the present case. On application of the 
field, the local magnetic moments become co-polarized 
giving rise to a FM state. The difference in the magnetic 
ordering can affect the conductivity and induce a mag- 
netoresistive (MR) response as discussed for graphene 
nanoribbons in Ref. llSl . To study the possible magne- 
toresistance signals in hydrogcnated graphene samples, 
the spin density and related spin-dependent on-site ener- 
gies have been computed for this specific case [Figl2jc,f)]. 

Kuho conductivity methodology - Using the self-consistent 
Hubbard Hamiltonian calculations, spin-dependent on- 
site energies are estimated for each supercell. These 
tight-binding parameters are then used for computing the 
conductivity in realistic-sized samples applying the Kubo 
formalism. These graphene samples are constructed by 
attaching a large number of these supercells following 
different random arrangements. An efficient real space 
order N method is then used to simultaneously follow 
the wavepackets dynamics and to compute the Kubo 
conductivity (see Refs. for details). The analysis of 
the wavepacket dynamics allows to extract the conduc- 
tion mechanisms and transport length scales (mean free 
path). As long as < 0.5% the diffusion coefficients are 
found to reach a saturation regime after a few thousands 
of femtoseconds, and then remain roughly constant up 
to very long times, indicating a diffusive regime and a 
negligible contribution of quantum interferences within 
the computational reach. The maximum value allows 
the evaluation of £e = Dmax/vp, with vp = lO^ms"^. 
Note that the transport properties are calculated in the 
zero temperature limit, which should however be robust 
to higher temperatures given the weak electron-phonon 
coupling in graphene. 

Spin- dependent transport features.- The scaling analysis 
of the conductivity is carried out using the Kubo for- 
mula cT^^iiE,t) = {ey2)Tv[S^^^{E - H)]D^,^{E,t) where 
Tr[^-|-4(£^ — H)/S] are the spin-dependent DoS per unit 
of surface at the energy E. At a fixed energy, tj-|-.4^(£^, i) 
reach their maximal values at the same time as D-f^i(E, t) 
(D™^^) which gives the spin-dependent Drude conduc- 
tivities aPl^'^'^iE) = (eV2)Tr[Jt,i(-S " H)]D'!['f'{E)), 
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FIG. 3. (color online) Left panel : Conductivities of hydro- 
genated graphene samples with various concentrations for the 
antiferromagnetic case (left panel), paramagnetic case (right 
panel inset) and ferromagnetic case (right panel main frame). 
Dashed lines give 2e^ /-kH. 



that accounts for the disorder effects on the density of 
states. 

Fig. [3] presents the Drude conductivities for the AF (left 
panel), FM cases (right panel) and the approximated FM 
case (inset in right panel) for H concentrations varying 
from 0.1% up to 0.5%. As expected, both in the AF and 
the PM case, no spin- dependent transport features are 
observed in the results, with the conductivities for up 
and down spins being equal within the statistical error. 
We note, however, some differences in the shape of the 
conductivity curves versus energy close to the Dirac 
point between the PM and AF cases. 

In sharp contrast, in the FM case, the spin-dependent 
conductivities significantly differ away from the charge 
neutrality point, with |cr°™'i'=(£:)/cr]^™'^<'(£;)| typically 
increasing with energy. a^^^'^°{E) decays roughly 
linearly with Ux for the lowest H density (which is a 
trend general for resonant scattering [2l|), but one 
clearly observes a tendency towards saturation at higher 
concentration. It is worth mentioning that, in all cases. 



the total Drude conductivities (cr|^™^'=(£') + cr|^™^<=(£')) 
remain larger than 2G'o/7r = Ae'^/irh, which is the 
theoretical minimum value estimated within the self- 
consistent Born approximation [2^ . 

^^From the computed conductivities in the AF 
and FM cases, one extracts (for AB-sublattice 
preserved symmetry), the magnetoresistance 
^^-^ "a^ + a^^), with = + af 
j-f-^. Fig. H shows MR for three 
different H concentrations. The maximum signal is 
observed for the lowest density Ux = 0.25% and reaches 
about 7% for a switch from the antiferromagnetic to 
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the excited ferromagnetic case. The signal is reduced 
with increasing and tends to disappear at higher 
concentrations, along with the spin density. One also 
notes that, since the macroscopic states are built out of 
supercells carrying only two vacancies, the number of 
magnetic correlations is likely underestimated. Hence, 
experimental observations of larger MR values would 
not be a surprise. We observe some asymmetry in 
the MR profiles (Fig. 4) which roots in small numerical 
inaccuracies stemming from both, the approximated 
Hubbard Hamiltonian and the use of a random phase 
approximation for the Kubo methodology. Fig. 4 (inset) 
also emphasizes the energy- dependence conductivities 
for Hx = 0.25% and = 0.35% close to Fermi level, 
which result in the qualitatively different profiles ob- 
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FIG. 4. (color online) Magnetoresistance of hydrogenated 
graphene samples with concentrations of = 0.25, 0.35, 0.5% 
as a function of charge energy. 



Conc/wsion.-Spin-dependent conductivities in disor- 
dered, hydrogenated graphene have been computed 
and related to the underlying magnetic ordering taking 
place at a large scale. Magnetoresistance signals of 
up to 7% have been found in the low density limit, 
suggesting that the existence of local magnetic states in 
graphene, induced by hydrogenation, could be observed 
experimentally. The high concentration limit [ux > 1%) 
in which stronger quantum interferences and localization 
effects should take place deserves a further consideration. 
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